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'. Abstract 

(N 

Generalized models provide a frame"work for the study of evolution equations with- 
^/^ . out specifying all functional forms. The generalized formulation of problems has been 

Q I shown to facilitate the analytical investigation of local dynamics and has been used 

successfully to ans"wer applied questions. Yet their potential to facilitate analytical 
computations has not been realized in the mathematical literature. In the present 
paper we introduce the method of generalized modeling in mathematical terms, sup- 
porting the key steps of the procedure by rigorous proofs. Further, we point out open 
questions that are in the scope of present mathematical research and, if answered could 
greatly increase the predictive power of generalized models. 
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1 Introduction 

Many processes observed in nature are too complex to be described on a detailed mechanis- 
^ , tic level. Therefore mathematical modeling can often not provide an exact set of evolution 
^ ! equations. This particularly evident in the context of mathematical biology j24] : 

"The Hodgkin-Huxley models are based upon sound biophysical principles, but these prin- 
ciples do not constrain the models to a definite set of equations [. . .]. " 

Therefore we always seem to face the dilemma that a dynamical analysis requires a given 
specific model. However, once we specify some of the functions, that are only partially or 
not at all known, then we cannot provide a result that is valid for the underlying physical, 
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chemical or biological process in full generality. 

One possibility to address this problem is to introduce a wide variety of parameters or 
even phase-space variables in an ad-hoc way. A more systematic approach is provided by 
the theory of S-systems [8] that aims at grouping different terms in evolution equations. 
Another systematic approach is considered in metabolic control theory [371 [3T| H6] where a 
linearized analysis for a dynamical system uses sensitivities as a standard set of parameters. 

Generalized modeling applies normalizing coordinate transformations to a system grouped 
into gain and loss terms to obtain a systematic parametrization. A generalized model pro- 
vides an intermediate alternative enabling the mathematical modeler to use the partial infor- 
mation he has available but still provides enough flexibility to treat many different alternative 
models simultaneously. 

Sections [2ll3] form an extended introduction to readers not familiar with generalized mod- 
eling whereas we focus on new results in Sections HHHl In Section [2] we review several results 
obtained in applications to show what conclusions can be drawn from a generalized model. 
For illustrating the application of the method in practice we briefly consider an example of 
a planar predator-prey model in Section [31 

In Section |4] we provide systematic treatment of generalized models for arbitrary ordi- 
nary differential equations (ODEs) in M". We analyze the normalizing transformations for 
generalized models in detail, supporting this key step of the procedure by mathematical 
proofs. We also prove results on the number of parameters and discuss the role of positivity 
assumptions. This approach provides a standard scheme for the application of generalized 
models under minimal mathematical assumptions. 

In Section |5] we answer several mathematical questions that arise from the framework of 
generalized modeling in the context of results and methods from bifurcation theory. The fo- 
cus of this analysis is on genericity and the structure of bifurcation diagrams. We explain via 
several instructive examples for non-degeneracy conditions what information can or cannot 
be inferred from a bifurcation analysis of a generalized model. 

In Section [6] we show that generalized models can also be applied to homogeneous steady 
state dynamics of a wide variety of other evolution equations including delay, partial and 
stochastic differential equations. We also indicate some recent progress on extending gener- 
alized modeling to nonlocal dynamics for the case of periodic orbits. 

2 Results in Applications 

Generalized modeling was originally proposed in the context of community ecology [21]. Only 
subsequently it was recognized as a general approach [22] and applied to a wide verity of 
apphcations. A comprehensive list of publications is[2Il[221l2SlHISlESlESll5Sl[ISl[ISlE51 
[T7t [Ml [511 |50l [511 [m [3l [131 EI] ■ In these apphcations generalized models for instance revealed 
essential factors for the stability of food webs [23] ; resolved a discrepancy between different 
modeling approaches in ecology [51]; implicated a Hopf bifurcation as a cause of Paget's 
disease in humans ^T\] and provided insights in the stability of mitochondrial metabolism 
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Besides the ability of generalized modeling to deal with unspecified relationships, the 
success of the approach builds mainly on making the Jacobian matrix analytically accessible. 
As shown in more detail below, the Jacobians obtained form generalized models are given in 
explicitly and typically contain only simple functions of the parameters. In systems of small 
and intermediate size the bifurcations of a generalized model can thus often be computed 
explicitly by hand. 

For exploring larger generalized models one typically randomly samples the local stability 
in random points of the parameter space to build up a database that is subsequently explored 
by machine learning techniques. Because of the direct accessibility of the Jacobian matrix 
the evaluation of every sample point involves only the computation of the leading eigenvalue 
of a matrix. Because of the numerical efficiency of this computation the exploration of large 
parameter spaces becomes feasible. For instance in [23] a system containing 50 dynamical 
variables and thousands of unknown parameters was explored by 100 billion (10^^) samples, 
which were obtained in reasonable numerical time. Generalized modeling is thus one of very 
few approaches that has reasonable hope of scaling for instance to whole-organism models 
in systems biology. 




Figure 1: Bifurcation diagram in generalized parameter space. Re-printed with permission 
from [52] ; see also |55]- The red surface indicates Hopf bifurcation and the blue surface 
saddle-node bifurcations. Higher co-dimension bifurcations are indicated by separate labels. 
The main codimension two curves are Gavrilov-Guckenheimer (GG), Takens-Bogdanov (TB) 
and double-Hopf (DH) bifurcations. 

A major limitation of generalized modeling - the restriction to dynamics close to equi- 
libria - is so far mitigated in applications, by the analysis of local bifurcations of higher 
codimension. These bifurcations allow some conclusions on global dynamical properties and 
can be computed analytically in generalized models up to intermediate size and numerically 
beyond that. Figure [T] shows an example of the information that is obtained. The three- 
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dimensional bifurcation diagram shows codimension-one bifurcation surfaces (fold and Hopf ) 
as well as codimension-two curves (Gavrilov-Guckenheimer, Takens-Bogdanov and double- 
Hopf) and codimension three points (1 : 1 resonance, triple point Although not all 

the unfoldings of the higher-codimension bifurcations are known it often suffices to detect 
the bifurcation point as an organizing center to gain insight into the overall dynamics. For 
example, parts of the double-Hopf bifurcation are known to generate chaotic dynamics due 
to associated torus and homoclinic bifurcations |25] . 



3 A Predator-Prey Example 

To illustrate the basic steps in generalized modeling we apply it in a non-rigorous way to a 
planar predator-prey system [22] with a prey density X and a predator density Y. The prey 
population grows at a rate S{X), predation occurs at rate G{X, Y) and natural mortality of 
the predator at rate M{Y) which yields 

X' = s(x)-G(x,y), 

Y' = G{X,Y) - M{Y). 

Suppose ([1]) admits an equilibrium point {X,Y) = {X*,Y*) and introduce normalizing 
coordinates 

X ^ Y 

X := — and y := —. (2) 

moving the equilibrium to {x,y) = (1, 1). Then we normalize the rate functions 

SiX*x) GiX*x,Y*,y) M{Y*y) 



Direct substitution of ([2])- ([3]) into ([T]) gives 

x' = ^^s{x) -^^^^^g{x,y), 
y' = ^^^y^r^ g{,x,y) - ^^^m{y), 

where we define new parameters 



Since (x, = (1,1) is an equilibrium point we know that the following holds: 

= ^^^7(l,l)-^m(l) = 
Therefore (HI) can be re-written as 



(4) 



_S{X*) _ G{X\Y*) _ G{X\Y*) _ Mix*) 

P-" ■- j^* ' ■- J^* ' •- ' P"" ■~ Y* ' ^ ' 



(6) 



x' = (3i{s{x) ~ g{x,y)), 
y' = (^2{9{x,y) - m{y)). 
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The Jacobian at the equihbrium {x,y) = (1, 1) is then given by 

V /32 5x[^(a;,?/)]|(:.,y)=(i,i) f32 dy[g{x,y) -m{y)]\^^^y)=i^i 

(^29x P2[9y-my] J 
where dx, dy denote partial derivatives and we have introduced another set of parameters 

Sx = dx{s{x))\x=l, gx = dx{g{x,y))\(^x,y)=il,l), qqn 

gy = dy{g{x,y))\(^x,y)={i,i), ruy = dy{m{y))\y=i. 

Based on this parametrization of the Jacobian, one can start to carry out a bifurcation 
analysis. From this example, we observe that the main steps of the method are: 

(Gl) Build a model of the underlying process as evolution equations and group different 
terms in the resulting equations. 

(G2) Apply a transformation in phase space based on the existence of an equilibrium point 
and introduce parameters. 

(G3) Interpret the generalized parameters in the modeling context. 

(G4) Apply methods such as bifurcation analysis to characterize the dynamics of the gener- 
alized model. 

However, many open mathematical questions remain when we try to apply (G1)-(G4). 
It is the main goal of this paper to provide a detailed systematic and rigorous description of 
(G1)-(G4) for a wide variety of evolution equations. Furthermore, we will answer many ques- 
tions that could not be answered appropriately in the context of applications; see references 
in Section [21 



4 Structure of Generalized Models 

We start with generalized models for ordinary differential equations (ODEs). A general 
autonomous first-order system of ODEs is given by 

— = X' = F{X;f^) (11) 

where X e M" are phase space variables, n E are parameters and the vector field 
F : X — )■ is assumed to be at least continuously differentiable in X and continuous 
in fi. If we can specify a particular map F then the main task is to analyze the dynamics of 
(ITT]) i.e. to partition the parameter space MP into regions of qualitatively equivalent dynamics 
[13]. If we do not specify any assumptions on F we focus on the abstract analysis of ODEs 
[28] . Generalized models provide one possibility to bridge the gap between specific models 
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and abstract analysis by making some structural assumptions on F without specifying the 
map completely. We assume that fllip has a decomposition of the form 

X[ = = J]a,,fcF,,,(X;/i), (12) 

k=l 

where ai^k = +1 or aj_fc = —1, the subscript i G {1,2, ... ,n} indicates the i-th coordinate, 
Ki > 2 and Fj, ^ : M" x — )■ M+. We are going to discuss the positivity assumption on 
the map F at the end of this section. The terms Fj ^ with Oj fc = 1 are called gain terms and 
those with Oj^fc = —1 loss terms. We note that the type of decomposition is decided as part 
of the mathematical modeling and does not follow a fixed set of rules. However, the basic 
principle of grouping the different terms is often provided by their role in the mathematical 
model and the resulting system (fT2l) has a systematic structure. 

A first step to understand the dynamics of f|TT]) is to analyze the stability and bifurcations 
of equilibria. Suppose there exists an equilibrium point X* so that F{X*; fi) = 0. If we are 
only interested in the local dynamics near X* we can relax the differentiability assumptions 
on F to a neighborhood of X*. The local dynamics at X* is given to first-order by analyzing 
the eigenvalues of the Jacobian 

^(^;-")lx=x* = PxF)(X*;/.)= . (13) 

If we do not specify F exactly then X* has to be treated as an unknown. The derivatives 
of functions/rates at the unknown equilibrium are often difficult to interpret in terms of 
physical parameters. Therefore, we would like to consider a transformation that allows a 
physical interpretation of parameters. Generalized modeling [21] assumes that X* 7^ 0; we 
are going to discuss the special cases X* = and X* — at the end of this section. Then 
one considers the normalizing coordinate change 

x^ = ^=■. h,{X), for 2 G {1, 2, ... , n}. (14) 

We remark that the idea of re-scaling to simplify or de-singularize a problem appears in 
several mathematical approaches for analyzing nonlinear systems; a typical example is pro- 
vided by the blow-up method [5l[6] that can be viewed as a phase space re-scaling in suitable 
coordinates. With the transformation (fT4l) the ODE (fTTl) transforms to 

x'i = ^FiiX*Xi,X;x2,...,X:Xn)=:Fi{x), for z G {1, 2, . . . , n} (15) 

where we have omitted the /i parameter dependence for notational convenience. Therefore 
the equilibrium X* is transformed to x* = (1, 1, . . . , 1) =: 1. 

Remark: A standard coordinate change in dynamical systems [25] is to consider the 
transformation Xj := Xj — X* so that the equilibrium point is moved to a; = (0, . . . , 0) =: 0. 
This transformation is mathematically convenient but does not provide a normalization of 
parameters as the coordinate change (fl31) . 

We can immediately check that the eigenvalues of the Jacobian 0131) remain unchanged. 
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Proposition 4.1. If X* ^ for all i & {1, ... , N} then the eigenvalues of the Jacobian are 
invariant under (IT^ i.e. spec{DxF{X*)) = spec{D {!)) . 



Proof. By direct calculation we find that 

d 



det 



(D,F)(l)-AId 



det 



det 



dxj XI 



F,(X*xi,X*a;2,...,XX) 



-A Id 



x=l 



X*dX, ' X* ^ 



(a) 



fl- 



det 



vj=l 
n 



(9F- 



\j=l / \i=l * 



det 



dX, 



= det[(Dxi^)(^*) - A Id] 

where we have factored out non-zero scalars in step (a) for each row and in (b) for each column 
using linearity of determinants with respect to rows and columns. The result follows. □ 

The equivalence of eigenvalues and the associated stability properties turns out to be 
of primary importance in many applications of generalized models [23]; see also Section |2l 
However, Proposition 14. II can also be viewed as a corollary to the following global result. 

Proposition 4.2. Suppose that F & for some A; G Nq U {w}, where denotes analytic 
functions, and that X* ^ 0. Then the ODEs (ITT!) and (fT5|) are C^-smoothly equivalent via 
the map 



Proof. Observe that x = h{X) = {hi{X), . . . , hn{X)) is a C'^-diffeomorphism that conjugates 
the vector fields 

F{X) = {Dxh)-\X) F{h{X)). 

□ 

The global smooth equivalence we showed is much stronger than (local) topological equiv- 
alence |l3] and the normalizing coordinate change (fT^ can be viewed as leaving the dynam- 
ics completely unchanged. The next steps of generalized modeling involves grouping and 
labeling the free parameters so that they can be interpreted as modeling parameters. We 
introduce a notation for the normalized gain and loss terms 



fiA^) 



Fi^ki^lXi., . . . , X^Xn) _ -Fjjc(X*Xi, . . . , X*x„) 



Fi,k{X* 



(16) 



i,k 



where 



Fi^k{X*) and we assume that Fi^k{X*) ^ 0. Then ([T5]) reads 
F*, 



k=l 



XI 



-fi,ki^), 



for z G {1,2, . . .,n}. 



;i7) 
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As a next step we group two parameters together 

kk ■■= §7- (18) 

Remark: From a physical point of view, we can also consider the units in the definition 
of Pi^k- F*^ is always a rate, for example mass per unit time. Since X* has the dimension of 
mass this implies that has the dimension 1/time. 

Using definition f[T^ and the equilibrium point condition = at x* = 1 in ([T7|) gives 
n conditions 

= ^ai,fcA,fc for i G {l,2,...,n}. (19) 

k=l 

Therefore, we can hope to eliminate n parameters. For example, we could try eliminating 
Pi^i and set Pi^i = aj^i(— Pi^k + X^a ^.=-i Pi,k)- This elimination procedure can be 

formalized as follows: Define the vector 

/3 := (A,i, /9i,2, . . . , Kk,. • • • , kKr^f e (20) 

where k = Yl^=i ^i- Then ( !T9|) can be re- written as a matrix equation 

Q = AP (21) 

where the n x k matrix A has elements in { — 1, 0, 1}. The rank-nullity theorem gives 

K = dim(ker(^)) + dim(im(^)) = dim(ker(^)) + rank(^)). (22) 

This shows that rank(^) is the number of parameters that we can eliminate and that 
dim(ker(^)) is the number of remaining parameters after the linear relations f ll9p have 
been applied. The elimination of parameters is related to concepts used in structural kinetic 
modeling where the matrix A is closely related to the stoichiometric reaction matrix with 
normalized entries |5Tl |32l EO] . 

Observe that we have a choice which parameters Pi^k we eliminate using the algebraic 
equations ( 12T|) . A further optional step is to introduce a parameter for each variable and 
set 

:= (23) 

where we assume that a, > 0. 

Remark: From a physical point of view, we want to introduce to nondimensionalize. 
This implies that has to have the dimension 1/time while (3i^k is dimensionless; in this case, 
we can interpret /3i^k as ratios. One particular important choice to make this interpretation 
more precise is to consider the possible definition [2T] 

«i:= kk= kk (24) 
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where the equahty between the two sums follows from (12T]) . Using this definition we find 
that 



i,k 



Xlfe: a, fe = l (^i,k Ylk: a, fe = -l 



which interprets (3i^k as the rate associated to the term with index {i, k) divided by the total 
gain (or loss) rate i.e. we have obtained a ratio; see also Section [2j We shall not make 
explicit use of definition ( 12^ here as it can be viewed as one particular choice of nondimen- 
sionalization. One can define a time scale as a physical quantity that has units 1/time. 
Therefore we shall call ctj time scale parameters from now on. Note that this justifies 
our assumption > on the basis of the underlying physical process. 

Now we can re-write the differential equation (fTTIl as 



Ki 



^ai^k/3i,kfi,k{^) ] 1 for i G {1,2, . . . ,n} (25) 
.fc=i / 

where the relation ( |2T]) is understood to apply as well. We call the parameters and Pi^k 
(resp. Pi^k) scale parameters. Obviously we have introduced quite a number of scale 
parameters to avoid specifying the functions in our model; therefore it is important to know 
how many scale parameters will appear in the model. We have the following result: 

Proposition 4.3. The number of scale parameters for a generalized model (|T71) is as follows: 

(CI) If Pi^k ore the only scale parameters and Pi^k 7^ A,m for all pairs {i, k) ^ (/, m) then the 
minimum number of scale parameters is given by 

K — rank{A) = dim{ker{A)) . (26) 

// all Pi^k appear as multiplicative factors after the elimination via Ap = then one 
more parameter can be eliminated. 

(C2) If Pi^k, Oii are the scale parameters and /3i^k 7^ A,m for all pairs (i, k) ^ {l,m) then the 
minimum number of scale parameters is 

dim{ker{A)) + n - T] - 1 (27) 



where rj is the number of scale parameters f3i^k that appear as multiplicative pre-factors after 
the elimination via f l?I]) . 

Remark: If the conditions Pi^k 7^ A,m resp. Pi^k 7^ Pi,m. are violated then further parame- 
ters can obviously be eliminated. However, a violation of this condition is not generic within 
the class of vector fields we consider here so we shall not consider this situation any further; 
for more on genericity see Section [51 
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Proof, (of Proposition I4.3p The previous discussion leading up to equation ( 122!) yields (126|) . 
The second part of (CI) that allows the elimination of one further parameter will be clear 
once we proved (C2). For (C2), we have dim(ker(^)) parameters (3i^k and n parameters a, 
after using the linear relations AP = 0. Assume without loss of generality that in the first 
r] coordinates, the parameters /3i^k, ■ ■ ■ , Pr),k appear as multiplicative prefactors so that the 
ODEs are 



"i/5i,iZ]fii«i,fcAfc(2;) 



(28) 



X' 



Now we define new time scale parameters ai := aif3i^k = Pi,k for i G {1, 2, . . . , r^} and a, = 
for i G {?7 + 1, . . . , n} which transforms (l28l) to 



x'l = ^lJ2k=l^^,kfi,k{x) 



x^ 



X 



'ri+l — Olri+lYlk=l f^ri+'^,kO'V+^,kfri+l,k{x), 



(29) 



n,k(^n,kfn,k \-^) ■ 



Therefore we have eliminated t] additional parameters. To eliminate one more parameter we 
can choose one time scale parameter a,, say without loss of generahty ai, and apply a time 
re- scaling 

t t/ai. 

Defining new parameters Oj/ai yields the final result. □ 

We continue by interpreting the parameters /Sj^fe. They are gain and loss ratios for each 
term in the decomposition. To analyze the stability of the equilibrium x* = 1 we define 
elasticities 

Ji,k,x, '■— — \ -L I 



dXj 



dxj 



[X 



(30) 



x=l 



The elasticities are sometimes also called or exponent parameters. Then we find the Jaco- 
bian of fl25l) at a; = 1 



J(l) 



/ ai 

^2 



\ 



V 







Or 



/ J2k=l^^,kf^l,kfl,k,xi ■■■ ^k=l^i,k/3l,kfl,k,x„ \ 



\ y^yfr— ] (^n,kPn,kfn,k,x\ ' ' ' 5^fc=l ^n,k(^n,kfn,k,Xn ) 



■ (31) 



We also refer to the set of scale and exponent parameters as generalized parameters. 
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The key input to the generahzed modehng process from apphcations is that we can often 
interpret the scale parameters and the elasticities for a given application. To explain the 
term elasticities we examine their definition more closely and observe that 



dx. 



^(1) 



lnF,,fc(X) 



x=x* 



which interprets the exponent parameters as (scaled) logarithmic derivatives. Logarithmic 
derivatives are often called elasticities, particularly in the context of modeling economic 
problems [TT] and in metabolic control theory [TU] . We can also view the exponent parameter 
fi,k,xj as the sensitivity to variations of fi^k in the direction Xj at the equilibrium point 
[l9]. If certain specific functional forms Fi^k{X) are known from the modeling process we 
can get even more information (see [16], p. 49). We give a few examples using uni-variate 
functions Fi^k{X) with X G M: 



F^,k{X) 


f 


AX'i 

Aexp(5X«) 

A 

B1+X1 
AXP 


qBX* 

(^B+{X*)P){Bp+ip-q){X')1) 


B+X1 





Further possible dependencies of the exponent parameters are easily derived by direct 
differentiation of the given functional form. 



We also have to make sure that we can find generalized parameters that belong to at 
least one specific model. Here, we support this important condition for sampling analysis 
[23]. 

Proposition 4.4. Suppose we are given scale parameters that satisfy ( !T9|) . Given any set 
of elasticities fi^k,xj there exists a specific model for the given generalized parameters. 

Proof. The proof is constructive. We are going to define functions Fj ^(X) that will produce 
the given set of generalized parameters. Fix some i G {1, 2, . . . , ra} and k G {1, . . . , Ki\. 
Then define 

where we can assume without loss of generality that 'ji^k 7^ and we can choose qi^kj for 
J G {1, 2, ... , n}. It follows from ([30]) that 



fiA^) 



%,k{x*y^.^'^xf'''''{x;y^.^'^xl''''-' ■ ■ ■ {x*y^^^'^x'^:'''-" 



7i,fe(Xr)«^M(X2*)''^fc.2...(X*)''«.^. 
Calculating the elasticities yields 



/■ k 



dfi,k ^-^-^ 
dxi 



1i,K 
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Since we are free to choose qi^kj we can match the prescribed elasticities. Note carefully that 
the previous calculation was independent on the choice of ji^k- Now define 



and observe that 



I?* f V* v*^ 

— Pi,k 



XI 

so that we also match the prescribed scale parameters. □ 

As a last step we re-visit the positivity assumptions on the maps F and the assumption 
X; ^ for the equilibrium X* = (X*,...,X*). First, suppose that Fi,fc(X*) = i^*^ = 
which implies that the definition ( fT6l) cannot be used to define the function /j^fc- In this case, 
we have to replace each term of the form Pi^kfi,k,xj in the Jacobian (13T1) by the standard 
un-normalized term 

However, suppose we know from the mathematical modeling that 

n 

F,,(X) = 5^(X,-XP''' 
1=1 



is a power function for some q = (gi,...,g„) with qi > 1 for all / G {!,..., n}. Then 
[f3i,kfi,k,x,] = l/c(i if qj = 1 and [Pi,kfi,k,xj] = if qj > 1. Hence we can discard the terms 
arising from Fi^k in the Jacobian if we know that ^ is locally super-linear and vanishes 
at the equilibrium. Note that the Jacobian is undefined for < g < 1 but that this case 
is also excluded by our assumption that F is at least continuously differentiable near the 
equihbrium point. 

Next, we consider the situation when X* = and assume that Fj_fc(X*) ^ for all i, k. 
Then the normalizing re-scaling transformation f|T^ is not well-defined. If we know from the 
mathematical modeling that X* = regardless of the parameters G M"*^ in the underlying 
model then we can simply define 

Xi . Xj 

and carry out the generalized analysis. If X*(yu) — )■ as /i — ?■ /i* for some fi* G MP then this 
situation can also be incorporated into the generalized analysis. Indeed, notice that 

n ^i,k M^/t* 1 r X* dFi k fJ.->-tJ* „ 

A,. = ^ ^ - = j-^^^ix ) ^ 0. 

Therefore, the cases of large scale parameters (3i^k (or /3j fe) and small elasticities involving 
partial derivatives with respect to Xi naturally incorporate the cases where X* tends to zero. 



12 



5 Bifurcations of Generalized Models 



In the previous section we have focused on the general algebraic structure of generalized 
models and the normalizing transformations. The next step is to investigate bifurcations in 
generalized models. We start with a brief review of some basic terminology and results that 
we are going to need throughout our analysis. 

Using the n x n Jacobian matrix fl3T]) we have access to the eigenvalues and their mul- 
tiplicities at the equilibrium point x = 1 {X = X*). One can calculate the eigenvalues 
Aj numerically using standard methods such as un-symmetric QR factorization [13] . If 
3ft(Aj) 7^ for all i G {l,2,...,n} the Hartman-Grobman Theorem [SU] implies that the 
flow near x = 1 is locally topologically equivalent to the flow of the linearized system. In 
particular, we get asymptotic stability if 3fJ(Aj) < for all i. Therefore a necessary condition 
for bifurcation under parameter variation is that 3fJ(Aj) = for one (or multiple) eigenvalues. 
We briefly recall how to define an unfolding of a generic bifurcating family [571 12] as we need 
this terminology throughout this section. Let F G C^(]R",M"') be a vector field defining the 
ODE (|TT1) . The smoothness r will not be of primary relevance for us and we always assume 
that r is sufficiently large in the following, at least locally. Let J^iF) denote the s-jet 

{x, F{x), DF{x), D'^F{x), . . . , D'F{x)) 

of F at X with s < r; denote the associated space of jets by J*(]R",M"). The s-jet extension 
F of F is a map 

F : R" J'(R", M"), F(x) = J^(F) 

that maps a phase space point to the associated jet; observe that we can identify the jet 
space J''(R"',]R") with R™ for a suitable m. Let Xq G R" denote a phase space point and 
let U he a neighborhood of Xq and set V := F{U); we shall restrict to studying the local 
behavior of the vector field near xq from now on. Let S C J'*(f/, V) denote the codimension 
n subset of those s-jets that have an equilibrium point xq in U [57J where codimension is 
defined as 

codim(^^) = dim(J*(f/, V)) — dim(£^) = m — dim(£^). 

Let B G S G J^{U, V) denote the codimension n + 1 set of vector fields with a non-hyperbolic 
equilibrium point xq; observe that a non- hyperbolic equilibrium point is defined by conditions 
on DF{xq). Consider an s-jet JJ(F) with a non-hyperbolic equilibrium point at xq. Then 
J^{F) lies in a set V G B of codimension b in J'*(f/, for b > n + 1 and we define 
the codimension of the equilibrium point Xq as 6 — n. To understand the dynamics near 
a bifurcation point in P we need a parametrized family of vector fields F(x; fi) so that 
the associated family of s-jets is transverse to V. Recall that transversality of maps G G 
C"'(R", R™) to a submanifold M G R*" at G{xo) is defined by the requirement 

F)G(xo)T,oR" + rG(.o)M = rG(..o)R™. 

The following theorem is of fundamental importance to justify the next steps. To state 
the theorem we recall that a property that holds for a countable intersection of dense open 
sets (i.e. on a residual set) is called generic. 
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Theorem 5.1 (Thorn's Transversahty Theorem, Let M be a submanifold of J^{U, V). 

The set of maps F G C'''{U, V) whose s-jet extensions are transversal to M is a residual set 
in C^{U, V) (for some r depending on s and n). 

Hence if we can find a parametrized family F{x] fi) with fi & MJ' which is transverse to 
P then we have constructed a generic representative. However, so far we have not taken the 
dynamics completely into account. Consider a vector field F* with F*{xo,i'o) = 0. We say 
that z/) is induced from near (xo,/io) if there is a continuous map 0, defined 

near fiQ with 0(z/o) = /^O; so that 

F*{x,u) = F{x,(f){u)). 

A parametrized family F{x,fi) is called a universal unfolding near an equilibrium point 
{xo,fiQ) if every other parametrized family of C" vector fields is equivalent to a family of 
vector fields induced by F{x,fi). In general, it is difficult to verify for many bifurcations 
with higher-dimensional parameter spaces that a transversal family also forms a universal 
unfolding. 

Remark: Although the results described so far give a framework for the classification of 
bifurcation points according to codimension there are a few subtle technical points regarding 
e.g. the applicability of Thom's Theorem 15.11 or possible re-parametrizations of time and 
coordinate changes |57]. 

The classification of local bifurcation according to codimension [13| [25] yields the follow- 
ing classification up to codimension two: 

• codim = 1: Fold or saddle-node (single zero eigenvalue), Hopf (pair of pure imaginary 
eigenvalues) . 

• codim = 2: Bogdanov-Takens (double real zero eigenvalues), Gavrilov-Guckenheimer 
or fold-Hopf (single zero eigenvalue and pure imaginary pair of eigenvalues) , Hopf-Hopf 
(two pairs of pure imaginary eigenvalues), cusp (single zero eigenvalue and degenerate 
normal form coefficient), Bautin or generalized Hopf (pair of pure imaginary eigenvalues 
and zero first Lyapunov coefficient). 

The eigenvalues of ( 13T|) depend on the generalized parameters so that we can detect 
necessary conditions for all codimension one and two bifurcations except cusp and Bautin 
bifurcation that are of codimension two as they violate codimension one non- degeneracy 
conditions. We shall not aim at a complete discussion of bifurcation analysis of generalized 
models but point out some features via examples. 

Example 5.2. Consider the following generalized toy model 

X' = Fi(X) - F2{X) + Fs{X) (32) 
for X G M. A possible specific model covered by f l32|) is 

X' = X^ -2X + 1- A, for AG (-1,1) (33) 
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where Fi{X) = X^, F2{X) = 2X and F3(X) = 1 - A. It is easily checked that ([33]) has 
two equihbria at X^ = 1 ± \/A for A > that undergo a non-degenerate fold bifurcation 
at A = 0. The bifurcation diagram is shown in Figure [2]^a). After the normalizing scaling 
transformation X = X^x we find that fl33l) gives two ODEs 

x' = Xlx^ -2x + =: 6±(x, A), for A e [0, 1). (34) 



X 



A 



1.1 r 

X 

1 

0.9 

1.1 

X 

1 

0.9 



1.1 



A 



(b- 



A 



0.9^ 
1.5 



2.5 



Figure 2: Comparison of conventional and generalized models in an example. Solid black 
curves are stable equilibria and dashed gray curves unstable equilibria, (a) Bifurcation 
diagram for the specific model (132|) showing a fold bifurcation, (b) Fold bifurcation for 
the re-scaled equations we plotted the two equilibrium curves at x = 1 separately to 
indicate that they are associated to different branches in (a), (c) Bifurcation diagram for 
the generalized model fl32|) . 

The new bifurcation diagram of flM|) is given in Figure [2](b±). The fold bifurcation from 
Figure |2]J^a) can still be recognized in Figure |2t^b±) at A = since 



dh 



1, 0) = and 



± , 



1,0) =2. 



dx ' dx'^ 
The transversality condition requires looking at the parameter derivative 



'dA 



[x,A) 



X 



2y/A 



(35) 



and we observe that the condition is not well-defined for {x^A) = (1,0). The generalized 
formulation of ( 132|) is 



x' = a{/3Mx) - ^2f2{x) + (/32 - l3i)Mx)). 



(36) 



Using a time re-scaling t i— )■ t/{a/3i) and setting (3 := /32//3i we get that (l36ll can be written 

as 



x' = /i(x)-/3/2(x) + (/3-l)/3(x). 



(37) 
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The associated Jacobian for ( 137|) at x = 1 is 

J(1) = /m-/3/2,x + (/3-1)/3,x. 

If we assume that the exponent parameters are (/i,x-/2,s, /s,^) = (2, 1, 0) as for our specific 
model then we find that 

J(l) = 2-/3 

which satisfies the necessary condition for a saddle-node bifurcation at /3 = 2. Since fl5B]) 
always has an equilibrium point at a; = 1 we find a bifurcation diagram as shown in Figure 
lie). 

From Example 15.21 we conclude that a generalized model will be able to detect the main 
bifurcation conditions for eigenvalue crossings but that one has to be careful in interpreting 
the shape of bifurcation diagrams. Furthermore, a first-order analysis will not check whether 
non- degeneracy and transversality conditions are satisfied. This raises the issue of genericity 
for generalized models which we address in the next example. 

Example 5.3. We continue with Example 15.21 and the generalized model 

x' = - PUx) + (/3 - l)/3(x) =: g{x,l3). (38) 

We check the generalized model for genericity near a fold bifurcation. Consider the associated 
1-jet space J^(]R, M) to g with elements 

x,g{x,l3),^{x,l3) 

which is three-dimensional. The bifurcation set B for a fold bifurcation is given by 

{(x,0,0) g4(M,M)} 

which has dimension 1 and codimension b = 2. Therefore we get, as expected, that the fold 
bifurcation has codimension 6 — n = 2 — 1 = 1. The 1-jet extension to (7 : R — > M is the map 

g{x,f3)= (^x,g{x,f3),^{x,f3)^ (39) 

which we view as a map from ^ : — t- including the generalized parameter /3. The 
linearization of the 1-jet extension is 



Dg{x,/3) 



/I \ 



(40) 



The tangent space to B is spanned by the vector (1, 0, 0)"^ G Tg(^x^p)M.^] note that we can 
obviously identify T^(^x,i3)^^ with M^. Evaluating ( HOl) at a fold bifurcation point (l,/3o) we 
know that 

|(1.A)=0. 
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This implies that checking transversahty of the 1-jet extension ( 139|) reduces to checking that 
the 3x3 matrix 



H := 

is non-singular. This just means 

det(if) 



1 








Oj 



1 \ 



0/ 



(9^ 



l,/5o) 



l,/3o)7^0. 



9/3 ' 

This recovers the well-known |l3] non- degeneracy conditions 

S'g 



and 



l,/3o)^0. 



(41) 



Recall that if fl4T|) holds then the unfolding is indeed universal in /3 [251 US] • Only the last 
step of the genericity analysis has to be adapted to the generalized model. In particular, we 
can plug in the structure of the model ( l38|) into ( 14TI) which yields the conditions 



|(l'/3o) 



l,/3o 



-/2(1) + /3(1)^0 



o&(l) + (/3o 



i)®(i)^o. 



(42) 



The conditions ( l42l) restrict the function space for which one can expect that the generalized 
model has a fold bifurcation with a universal unfolding. The examples have also demon- 
strated that known results about bifurcations and unfoldings carry over easily to generalized 
models. The same calculations can be carried out for all other local bifurcations which will 
yield conditions analogous to ( H2|) . 



Even if the unfolding is not universal (e.g. if we had chosen a function space in Example 
(15. 3 p with /2(1) = /3(1)) then we can often use the availability of additional parameters as 
the next example illustrates. 



Example 5.4. Consider a planar generalized model 



Efii«i,fcA,fc/i,fc(a;), 

Ef=i«2,fc/32,fc/2,fc(a;), 



(43) 



where we assume that the time scales are all equal to 1 for simplicity. Suppose fl43|) has a 
Hopf bifurcation at {x,Pi^i) = (l,/3*^) for all other generalized parameters fixed. The Jaco- 
bian J(l) at i = has a complex conjugate pair of eigenvalues Ai^2 = ^i,2{(3i,k, fi,k,xj) 
depending on the generalized parameters with ci := Trace(J(l)) = 0. There are two non- 
degeneracy conditions [13]. The transversahty condition is 



C2 



d 



dp 



-3?(A 



1,2, 



1,1 



7^0. 



(44) 



/3i,i=/3i*.i 



The second non- degeneracy condition is that the first Lyapunov coefficient li =: C3 is not 
equal to zero [25l HI] . The condition can be written in terms of the partial derivatives with 
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respect to Xi^2 up to third order. In total, we have to satisfy three algebraic conditions 
depending upon the generalized parameters and the choice of functions fi^k to get a non- 
degenerate Hopf bifurcation 

ci{(3i,k, fi,k,x,) = 0, (45) 
C2W^,k, h,k,.,) ^ 0, (46) 

C3i(3^,kJ^,k) ^ 0. (47) 

Observe that the scale parameters Pi^k or /32,fc will appear as coefficients in the linear com- 
bination of C3 so that varying one scale parameter we generically satisfy (H7|) in the space of 
smooth functions fi^k ^ C"'(M^,]R). Hence additional free generalized parameters can com- 
pensate for a restricted choice of functions fi^k- Varying two more generalized parameters 
we can always satisfy P5l) - fH6l) generically. 



We have seen that the scale parameters are relatively easy to understand. The elasticities 
are bit more complicated as the next example illustrates. 

Example 5.5. In principle, we can also just vary the elasticities and treat them as bifurcation 
parameters. The obvious caveat is that this might not yield a smooth family of functions 
everywhere in phase space. For example, if 

Fi,fc(X) = A{X - Bf, for p > (48) 

where A, B are parameters, then Fi^k is C°° for X ^ B but only for some finite /c if p ^ N. 
This can again lead to violations of non-degeneracy conditions for bifurcations as in Example 
(15. 2p . In particular, as shown in [IQ], we will not be able to conclude quantitative universal 
scaling laws near bifurcations if we do not restrict the functions ^ to be sufficiently smooth. 
Note that this phenomenon is again non- generic for a sufficiently large parameter space. The 
second key observation for elasticities relates to the form of /j fc. For fHHl) we find 

fi,k,x = _ (49) 

which depends on the equilibrium point location X* . The best way to interpret (1491) is 
asymptotic knowledge about X*; for example, we have 

X*>1, B = 0{1) h,k,.-p, 

X*^B, X* = 0{1) 

0<X*<1, B = 0{1) ^ Am~o. 

where 0(1) indicates a fixed constant independent of any asymptotic limits of X*. 

The last example shows that it can be beneficial to restrict the class of functions consid- 
ered in a generalized model to a certain class, e.g. smooth functions, polynomials or certain 
functionals particular to the application area [16j. 

The previous three examples illustrate the main aspects of the genericity question for 
generalized models: 
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1. We can check for a given decomposition and given function spaces whether the non- 
degeneracy conditions for a bifurcation are satisfied. 

2. The checking depends crucially on the underlying mathematical modeling and how the 
functions and decomposition are chosen. 

3. Higher-order terms have to be taken into account that are not parametrized by scale 
parameters and elasticities; see also [5U] . 

4. Bifurcation and non-degeneracy conditions can usually be satisfied upon considering 
sufficiently large generahzed parameter spaces. 

When studying an uncertain system one can not necessary conclude that a bifurcation 
detected by generalized modeling is crossed transversally when a given (conventional) pa- 
rameter is changed in the application. However, one can generically expect that bifurcations 
are transversal within the whole class of systems described by the generalized model. 

Another issue that was already briefly considered in (15.21) is the correspondence between 
generalized and specific models. If a specific system exhibits multistability one may ask how 
this is refiected in a generalized model. 

Example 5.6. Consider the specific bi-stable system 

X' = Q-11X + QX^-X^ (50) 

which is easily seen to have three equilibria X* = 1,2,3 where X* = 1,3 are stable and 
X* = 2 is unstable. A generalized model encompassing (150|) as a special case is 

X' = -Fi(X) + F2{X) - Fs{X). (51) 

We can obtain the special case by Fi(X) = IIX - 6, F2(X) = GX^ and F3(X) = X^. 
Also the generalized parameter sets representing (1501) can easily be found 

if X* = 1 (/3i, 132, (3,, /i,2,., /i,3,x) = (5, 6, 1, 11/5, 2, 3) 

if X* = 2 (/3i, /32, /33, fix-. /i,2,x, /i,3,x) = (8, 12, 4, 11/8, 2, 3) (52) 

if X* = 3 (A, /32, /33, /i,2,x., /i,3,x) = (9, 18, 9, 11/9, 2, 3) 

where we assumed without loss of generality for the single time scale that a = 1. The results 
fl52|) illustrate that three different points in generalized parameter space represent the three 
local equilibria. Indeed, we have considered generalized modeling here as a purely local 
method (but see also Section 16. 4p which also explains this representation of multistability. 

The previous example also raises a more general question about the correspondence 
between bifurcation results for specific and generalized models. In particular, one may want 
to relate the parameters /i G in ffTIl) to the generalized parameters which can be expressed 
via a mapping 

{/i G W} {ai, Pi^k, fi,k,x,}i,j,k- (53) 

Obviously this map M depends on the specific functional forms of the vector field F{x; /i). 
Therefore we shall not discuss the parameter space mapping here. Figure O illustrates the 
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Figure 3: Example of the correspondence between specific and generalized bifurcation dia- 
grams; figure adapted from [51]. (a) Bifurcation diagram for a specific model where /ii is 
a parameter and X2 a phase space variable. Thick lines denote equilibria with solid curves 
for stable and dashed curves for unstable equilibria. The thin lines mark the maximum and 
minimum X2 coordinates of a limit cycle generated via the supercritical Hopf bifurcation H . 
Two transcritical bifurcation TCi^2 and two saddle-node bifurcations are indicated as 
well, (b) Generalized parameter space bifurcation diagram where the two surfaces indicate 
Hopf (red) and saddle-node (blue) bifurcations. The curve represents the diagram from (a) 
obtained via a mapping (l53l) . To calculate the curve one follows the same procedure as in 
Example 15.61 

mapping M for a model system from ecology [51]. 

A very important conclusion of all previous examples is that generalized models produce 
information that "scales with the input" . The more knowledge from the modeling process is 
available, the more detailed information about the dynamics can be derived. 

We end this section by mentioning an important practical issue. We have to address 
how to find bifurcations and their associated varieties (such as surfaces, lines, curves, etc.) 
in parameter space. A direct method calculates the eigenvalues and uses iteration to find 
the zeros of real parts. An indirect method employs so-called test functions that vanish 
once a certain type of eigenvalue crossing occurs; see [15] or [13] for an overview of different 
test functions and references to the literature. We note that analytical methods such as 
the method of resultants [271 EH] have also been employed successfully in the context of 
generalized modeling [21] using computer algebra [20]; having an explicit formula for the 
bifurcation loci can be beneficial for visualization [531. 
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6 Beyond ODEs 



Although we presented the ideas of generahzed modehng only in the context of ODEs so far, 
it is evident that their scope is much broader. Recall that the key ideas are: 

• There exist unknown functional forms in the model. 

• Group the different parts of the vector field into gain and loss terms. 

• Introduce a normalizing coordinate change for an equilibrium point. 

• Re-scale the gain and loss terms and introduce scale parameters. 

• Linearize at the new equilibrium x = 1 and introduce elasticities. 

• Interpret the generalized parameters and assign suitable ranges. 

• Use tools from bifurcation analysis to capture the dynamics. 

These ideas carry over naturally into a much wider setup than just ODEs. Generalized 
models for discrete dynamical systems (iterated maps) are discussed in [39]. Many other 
mathematical evolution equations are very similar to ODEs. In particular, the notion of 
equilibrium point, coordinate changes or re-scaling as well as linearized analysis carry over. 
In this section we shall demonstrate this observation for certain classes of partial, functional 
and stochastic differential equations in Sections 16. 1116. 31 In Section 16.41 we provide a brief 
overview on extensions to nonlocal dynamics. 



6.1 Partial Differential Equations 

Generalized modelling for partial differential equations has first been developed for reaction- 
diffusion equations in the context of spatial predator-prey models [171 S] . A more general 
class of partial differential equations (PDEs) that can be analyzed by dynamical systems 
techniques is given by (semilinear) parabolic equations. Let {x, t) E Q x M+ where f2 C 
is a domain and let u = u{-,t) G X for all t > where X is a suitable Banach space e.g. 
X = LP{fl)] for details see [33]. Consider the abstract evolution equation 

Ou 

— + Au = F{x,t,u; ^) (54) 
at 

where F : x R+ x X x M*' — )■ R" and A is a sectorial (differential) operator on X [33]. A 
concrete example for fl54p are initial boundary value problems for reaction-diffusion PDEs 
[36] of the form 

^-VAu = F{x,t,u;fj.) forxefi, t>0 

u{x,t) = g{x,t) for X e t > (55) 

u{x, 0) = h{x) for X G dQ 

for g : Q X ]R+ — )■ M", /i : — )• R", A is the Laplacian and D is a diagonal matrix 
called the diffusion matrix with non-negative elements. We shall restrict our presentation of 
generalized models here to (|55l) but remark that the ideas naturally extend to many other 
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equations within the class ( 15^ . Assume that F can again be decomposed into gain and loss 
terms 

Ki 

Fi{u,x,t) = ^ai^kFi,kiu,x,t). 
k=i 

where we omit the parameters fi E W again for notational convenience. Note that this 
assumption is very natural in the context of reaction-diffusion systems since F represents 
the reaction terms. For example, if we think of a chemical reaction then gain terms would be 
terms that increase a concentration while loss terms decrease it. Suppose (j55|l there exists a 
constant function u* such that 

= K,...,0, , u*eR, <^0, 

for all t > 0. Then u* is a space-time homogeneous equilibrium solution to fl55l) . We can 
apply the normalizing coordinate change 

Vi := for i e {1,2,.. .,n}. 

Then the reaction-diffusion system ( l55ll transforms to 

Ef=i ai,k(3i,kfi,kix, t,v), xen,t>o 

X edQ,t>0 (56) 
xedn 

where the usual generalized modeling definitions as in fl25|) are used and the diffusion matrix 
has been rescaled Va = Va/ai. Then f l56|) has a space-time homogeneous equilibrium at 
V = (1,1,...,!) =: 1. Now all linearization techniques at f = 1 for flS^ can make use of 
the interpretation of scale parameters and elasticities as before. As a typical example one 
can consider the Turing- Hopf mechanism [36]. Let J(l) denote the Jacobian as in ( 13T|) . If 
J(l) has only eigenvalues A = \{k) with 3fJ(A) < then the equilibrium f = 1 is stable as 
a solution to v' = J{l)v. The Turing-Hopf bifurcation requires that some eigenvalue of the 
Laplacian has ^{X{k)) > for some k > 1 which results in a spatial instability and associated 
pattern formation. Explicit conditions for this scenario to occur can then be derived in terms 
of the generalized parameters and the rescaled diffusion coefficients [1]. 

6.2 Delay Differential Equations 

Generalized models for delay equations have first been considered in the context of coupled 
oscillators in |3l]. We are going to consider a slightly more general delay differential equation 
(DDE) with constant delays t/, I E {1,2, ... , m}, given by 

X' = F{X{t),X{t - n),X{t -T2),...,X{t- T^)- /i) (57) 

where F : x M?' -> M" and we abbreviate X{t - n) as X^'. As for PDEs, we remark 

that ( 157|) only presents a subclass of DDEs as it is non-neutral with constant delays ^29j and 



V,{x,t) = =] gl{x,t), 

v,{x,0) = !^=:h:{x). 
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we expect generalized modeling to apply for many more DDEs than just ( 157|) . Suppose there 
exists an equilibrium point X = X* so that 



F{X*,X*,...,X*;fi) = 0. 
Applying the usual generalized modeling procedure we end up with 



= 7 ,o.i,kl3i,kfi,kix) (5^ 



where the scale parameters are defined as usual and 

F,^k{X* ■x,X* ■x-\...,X* -x' 



X* 

The linearized equation for (|58|) around x = 1 is 



v' = J2 ^i^"' (59) 

1=0 

where the n x n matrices Ai consists of the rows 

y^Q»,fcA,fcy^(A/i,fc)(l) 



,fc=l 1=0 



and where Di denotes the total derivative of fi^k with respect to the l-th argument; note 
that we have employed the convention tq = so that / = denotes the first argument. The 
characteristic equation associated to (!59l) is obtained by assuming an exponential solution 
of the form e'^* and is given by 

det ^Ald - j = 0. (60) 

In contrast to ODEs, we see that fIBU]) is a transcendental equation which can have an infinite 
number of solutions A. It is known that if 3?(A) < for every solution then the solution is 
stable. Bifurcation analysis for (!60l) can then be carried out using numerical or analytical 
[9] methods. The generalized parameters are used as bifurcation parameters in this context. 

6.3 Stochastic Differential Equations 

Consider a system of stochastic differential equations (SDEs) [15] 

dX = F{X;ij)dt + G{X;iJ,)dW (61) 

where W = W{t) = {Wi{t),W2{t), . . . ,Wk{t))'^ is a fc-dimensional Brownian motion, F : 
M" — )■ R", G : R" — )■ M"^'^ is a matrix- valued function. The normalizing coordinate change 
can be applied to (1^ for a deterministic equilibrium F{X*) = with X* ^ for all 
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2 G {1, . . . , n}. Since the coordinate change is hnear the Ito formula [121 EH] reduces to the 
standard chain rule and we get 

dxi = • - ' XnXn, /i) ^^ _^ _l_(7.(Xia;i, . . . , for i G {1, 2, . . . , n} 

where Gi indicates the i-th row of G. The scale and exponent parameters for the deterministic 
drift terms Fi{X; fi)/X* can be defined as in Section |H We can also consider a normalized 
matrix- valued function for the diffusion term where each row is given by 

/ N Gi{XiXi, . . . , XnXn] IJ^) 

9i(x:N ■■= Y* ■ 

Then we can still formally linearize (I^Tj) at x = 1 and obtain to lowest order 



dxj 
^ k=i 



^dt + g{l; ij)dW (62) 

where ^ = C,{t) now solves an SDE with linear drift term and constant diffusion. There are 
multiple possibilities on how to develop a "stochastic bifurcation theory" [35j, \1\ . Therefore 
we shall not discuss generalized modeling for SDEs in any more detail. However, we expect 
generalized models for SDEs to work for bifurcation and stability analysis as well. 



6.4 Nonlocal Generalized Models 

So far, we have presented a systematic approach to generalized models to analyze local dy- 
namics. It is a natural question to ask whether the approach can be used to analyze invariant 
sets beyond equilibria such as periodic orbits, homoclinic orbits or tori. 

A first important point to notice in this regard is that generalized models already provide 
a way to locate certain global orbits via local bifurcations. For example, it is well-known 
[13] . under certain assumptions on normal form coefficients, that 

• Codimension one Hopf bifurcations imply the existence of a periodic orbit, 

• Codimension two Bogdanov-Takens points imply the existence of a homoclinic orbit, 

• Codimension two Hopf- Hopf bifurcations imply the existence of an invariant torus. 

In fact, under suitable conditions on equilibrium point for a saddle- focus homoclinic orbit 
or via the break-up of tori near Hopf-Hopf bifurcation one can infer the existence of chaotic 
invariant sets which are obviously global phenomena. This approach has been used for gen- 
eralized models in several applications [HI [22l [52], [55l |59l [61] ; see also Figure [H 

However, this approach cannot capture global bifurcations such as saddle-nodes of peri- 
odic orbits. Recent work of Kuehn and Gross [42j shows how to extend generalized modeling 
to periodic orbits for the predator-prey system ([T|). The main problem is that the general- 
ized parameter and fi,k,xj become time-dependent functions Pi^kif) and fi,k,xj(t). This 
changes the algebraic structure of the generalized model and introduces a so-called moduli 
flow constraint on the scale functions Pi,k{t)- The detailed description of this approach is 
beyond the scope of this paper and we refer the interested reader to [42j . 
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